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1. Introduction 

The properties of charmonium in the quark-gluon plasma is an issue of long-standing interest, 
following the suggestion [|l|] that suppression could be a probe of deconfinement. Lattice sim- 
ulations in the quenched approximation ^, ^ ^ have cast doubt over this picture, indicating that 
J /\\f may survive up to temperatures as high as 1.5 — 27^. Recently, potential model calculations 
using the internal energy of the heavy-quark pair have reached the same conclusion, and indicate a 
qualitatively similar picture in the case of Nf = 2 QCD 0]. This conclusion has, however, been 
disputed by other potential model calculations 

On the experimental side, data from SPS and RHIC show similar amounts of J/xjf suppression 
at both experiments, despite the big difference in energy density. Two different scenarios have been 
developed to explain this result. The sequential suppression scenario |^] takes its cue from lattice 
results, suggesting that the entire observed suppression originates from feed-down from the excited 
IP and 2S states, which melt shortly above T^, while the IS state survives up to energy densities 
higher than those reached in current experiments. On the other hand, the regeneration scenario 



[ ]10| , |11| , |12| ] suggests that additional charmonium is produced at RHIC energies from coalescence 
of c and c quarks originating from different pairs. 

Lattice simulations with dynamical fermions (2 or 2-1-1 light flavours) will be one of the es- 
sential ingredients in resolving several of these issues. In order to obtain a sufficient number of 
points in the euclidean time direction with moderate computational resources it is necessary to use 
anisotropic lattices, where the temporal lattice spacing a^- is much smaller than the spatial lattice 
spacing a^. Here we present results from Nf = 2 simulations with anisotropy ^ = as/a-j: = 6. 

2. Methods 

We have carried out simulations on lattices with a temporal lattice spacing « (7GeV)^' and 
an anisotropy ^ = a^/a^ = 6, giving a spatial lattice spacing of = 0.167fm. Configurations have 
been generated at a range of temperatures, from 220 GeV to 440 GeV, and two spatial volumes, 
8^(1. 34fm)^ and 12^(2fm)^. 500-1000 configurations were generated for most temperatures. The 
pseudocritical temperature for these lattices was estimated to be Tc =205-210 MeV. For details 



about the anisotropic lattice action, see []13[], and for further details about the high-temperature 



configurations, see | |14| ] 



The spectral functions Pr(ft>) are related to the imaginary-time correlator Gr(T) (we consider 
here only zero momentum) according to 

r dco , , , , , , cosh[a)(T- l/2r)l 

where the subscript F corresponds to the different quantum numbers. Correlators were com- 
puted in the pseudoscalar, vector, axial-vector and scalar channels, using local operators only, ie. 
r = {75,7i,757i,f}- Spectral functions were then determined using the maximum entropy method 
(MEM) JT5| , p^ ], with the modified kernel introduced in [p7||, which greatly improved the stability 
and convergence of the MEM analysis. 

In addition to the spectral functions, we have also studied reconstructed correlators Gied'^', T, Tq), 
obtained by integrating (|2.1[) at temperature T using the spectral function p{(o;Tq) obtained at some 



2 



Charmonium in the QGP 



J.I. Skullerud 




6 8 10 2 4 

(0 (GeV) (0 (GeV) 



Figure 1: The pseudoscalar (left) and axial- vector (right) spectral functions for different default models, at 
r = 221MeV. 

reference temperature Tq. If the spectral function at Tq can be obtained with a high degree of con- 
fidence, this may provide a more robust way of analysing possible medium modifications than by 
attempting to extract p{(o) using MEM at higher temperatures, where the number of data points 
decreases as \/T . 

3. Results 

The pseudoscalar and axial- vector spectral functions at T = 221 MeV, just above the decon- 
finement transition, are shown in figure i for a range of different default models. In the pseu- 
doscalar channel, we can see very little dependence on the default model: only for the default 
model m{(o) = 80g)^ is there any significant deviation. Note that in the graph, this model would 
corresponds to a horizontal line at 80, and in the absence of information from the correlators the 
MEM would reproduce this default model. It is natural that the tension between this rather ex- 
treme model and the data will produce some distortion in the reproduced spectral function. We 
conclude from this that the pseudoscalar spectral function can be reliably reproduced at this tem- 
perature. In particular, we note that the position of the first peak corresponds very closely to the 
zero-temperature r]c mass (the charm quark mass used in this analysis was somewhat lower than 
the physical charm quark mass). The second peak is assumed to be a lattice artefact as its position 
is close to a sharp cusp in the free-field spectral function. 

In the axial-vector channel, on the other hand, we find a great deal more variation. Our 'best' 
estimates for the spectral function (thick lines) indicate a surviving peak structure approximately 
at the mass of the zero-temperature Xci, but the strong default model dependence, combined with 
the proximity of the lattice artefact peak, makes it impossible to distinguish with any confidence 
between a bound state peak and a two-particle threshold. 

As discussed above, we can be quite confident about our pseudoscalar spectral function at 
r = 221 MeV, and will therefore use this as a reference temperature for studying reconstructed 
correlators. The same is the case in the vector channel, while in the P-wave (scalar and axial-vector) 
channels the spectral functions at T = 221 MeV are much more uncertain. We will nonetheless use 
this as our reference temperature also in this case. 
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Figure 2: Reconstructed correlators for the S-wave (left) and P-wave(right) channels. 
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Figure 3: Spectral functions in the vector (left) and axial-vector (right) channel, for different temperatures. 



The reconstructed correlators are shown in figure ^. We can see that in the pseudoscalar 
channel, Grec(T) changes by only 2% as the temperature is doubled with respect to Tq. In the vector 
channel the change is slightly larger, but still less than 10%. Although it is quite possible that a 
null result for Grec(i') can go hand in hand with dramatic changes in the spectral functions [^], this 
provides prima facie evidence that medium modifications are small in these channels. 

In the P-wave channels the picture is dramatically different: the long-distance correlators 
change by a factor 1.5-2 over a much smaller temperature range. This would indicate strong 
medium modifications, although a large part of the change is most likely due to the constant mode 



which appears in these correlators above Tc []18|]. This will be investigated further in the future. It 
is worth noting that our results for the reconstructed correlators are very similar to those obtained 
in the quenched approximation for the same values of T /Tc. 

Figure ^ shows spectral functions in the vector and axial- vector chamiel at a range of different 
temperatures. In the vector channel, the ground state peak appears to survive unchanged up to 
T 300 MeV, while atT ^ 440 MeV there is no sign of any surviving ground state. The default 
model dependence of our spectral functions becomes progressively stronger as the temperature is 
increased, and as a result we are not in a position at present to determine precisely where in the 
range T = 300 - 440 MeV theJ/xj/ melts. 

The nonzero intercept at ft) = at our highest temperature could be related to a non- vanishing 
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charm diffusion coefficient. We have attempted to study this using the default model m{(o) = 
moO}{b + (o) which explicitly allows limco^o p {(o) / (0 7^ [17]. However, our results at present are 
inconclusive, and finer lattices will in all likelihood be required to make any clear statements. 

The spectral function in the pseudoscalar channel is qualitatively similar to the vector channel, 
but the ground state appears to survive to somewhat higher temperatures. 

In the axial-vector channel there is no sign of any surviving ground state peak for T > 250 
MeV or 1.27^. Below this temperature there appears to be a peak structure which may correspond 
to a surviving Xci state or to a threshold enhancement as discussed above. We conclude from this 
that the Xci melts below 250 MeV, although a more careful analysis taking into account the constant 
contribution in the correlator will be required to draw any firm conclusion. A further complication 
is that the local operator has a poor overlap with the Xci state, making it more difficult to extract 
information about the ground state if it exists. This could be improved on by using a derivative 
operator in the axial-vector as well as in the scalar channel. The picture in the scalar channel is 
similar to that in the axial- vector channel, but the uncertainty is greater. 



4. Outlook 

Our results indicate that charmonium S-waves (J/YjTIc) survive in the quark-gluon plasma 
up to 1.4-2 times the critical temperature, while P-waves (Xc) melt below 1.27^. This is conis- 
tent with what has previously been found in the quenched approximation. There aie a number of 
uncertainties that must be addressed before firm conclusions can be drawn. 

• First and foremost, the spatial lattice spacing used in this first study is fairly large (0.17fm). 
This means that the number of temporal points in the relevant temperature range is quite 
small, even with an anisotropy of 6, limiting the usefulness and reliability of the maximum 
entropy method. Also due to the coarse lattice, lattice artefacts appear in the spectral func- 
tion uncomfortably close to the physical ground state masses, obscuring the physical inter- 
pretation of the spectral functions extracted using MEM. Simulations with finer lattices are 
underway, which should go a long way towards addressing these issues. 

• The simulation has used two degenerate flavours of sea quarks with a mass around the strange 
quark mass, rather than two light quarks and one strange quark as in nature. The main effects 
of this are to raise the (pseudo-)critical temperature above its physical value, and to inhibit 
charmonium dissociation by light quarks. Simulations with lighter sea quarks are planned. 

• The P-waves should be probed with operators which have a better overlap with the states of 
interest. We are currently producing P-wave correlators using derivative operators, which we 
expect will provide more robust results for these channels. 

We are also computing charmonium correlators at nonzero momentum. This will allow a 
more direct comparison with experimental results, where all J/y particles that are observed (as 
dileptons) in the detectors have nonzero (transverse) momentum. Most importantly, it may inform 
the interpretation of the data in terms of models such as the sequential suppression model and the 
regeneration model, since the momentum and rapidity dependence of the J/y yields is a key factor 
in distinguishing between these different models. 
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